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(D ' Abstract 



We consider scalar field theory in a changing background field. As an example we study the 
simple case of a spatially varying mass for which we construct the semiclassical approxima- 
tion to the propagator. The semiclassical dispersion relation is obtained by consideration of 
spectral integrals and agrees with the WKB result. Further we find that, as a consequence 
of localization, the semiclassical approximation necessarily contains quantum correlations in 
momentum space. 



1 Introduction 

The semiclassical (WKB) method has been developed to describe motion of a quantum 
particle in a slowly varying background. The method is a systematic expansion of the wave 
function in powers of the Planck constant, or equivalently, in gradients of the background [|l|. 
In field theory, on the other hand, one is usually concerned with coupling constant expansions, 
and applicability of the semiclassical method has appeared more limited. Nevertheless, there 
are many cases where the semiclassical method can be useful also in field theory. For example, 
it may be adequate to describe a part of the system by a classical field to which some other 
quantum field is weakly coupled. An important physical example of this occurs at the 
electroweak phase transition, where the problem is to study a particle propagation in the 
presence of a spatially varying, CP-violating Higgs condensate. Interactions of the plasma in 
this background are known to lead to baryon production via the electroweak anomaly. This 
is widely thought to be the most attractive explanation of the matter-antimatter asymmetry 
of the Universe [@], soon to be tested by accelerator experiments. 

The scale at which the Higgs condensate varies is determined by the coupling of the 
Higgs field to other species in the plasma and by the dynamics of the phase boundary during 
the phase transition P]. Quantitative studies have shown that phase boundaries are often 
thick when compared with a typical De Broglie wave length of a particle in the plasma [0. 
Because the Higgs condensate endows particles with mass, it is then reasonable to model 
the condensate by a spatially varying mass. It is known that source terms in the dynamical 
equation leading to baryogenesis, due to these mass terms, appear beyond leading order in 
gradient expansion |^, |^ . However, neither particle propagation, nor plasma dynamics has 
so far been studied in a controlled approximation scheme beyond leading order in gradients. 
Other applications of the semiclassical method abound, e.g. motion of electrons in the 
background of spatially varying potentials, etc. 

In this letter we develop an approximation in powers of gradients for the propagator of 
a scalar field with a spatially varying mass term m? = m?{t,x) in a theory described by the 
lagrangian 

C = (9^0)^(9» - m20t0 + £.^^ , (1) 



where £int contains interactions. This is an important ingredient required in a generahzed 
description of the plasma dynamics in a spatially varying background. To interpret our solu- 
tions we study the properties of spectral integrals of test functions for dynamical quantities. 
We argue that localization in space implies that the usual quasiparticle picture has to be 
extended to include a limited amount of correlations (quantum coherence information) in 
momentum space. 

Before considering the semiclassical method for the propagator, we derive the WKB 
dispersion relation from the wave equation. For simplicity we consider a stationary case, 
such that rn? = rn?{x) with x defined by x^ = (0,x, 0, 0). Ignoring £int, Eq. (|I]) then implies 

[uj^-kf,+dl-m\x))^^,:^^=Q, (2) 

where uj is the conserved energy and fc|| is the conserved momentum in the j/z-direction. 
Rewriting the wave function as $^^ = uexpi !^ kx{uj,k^(,x')dx' , Eq. (|]) can be recast as 
an equation for the momentum variable 

Here we used u' /u = —k'^/2kx and k^ = u'^ — k^ — rn? and e.g. k'^ = dxk^. Eq. (|) can be 

solved iteratively in the derivative expansion; neglecting for simplicity d^m?' for n > 2, one 

obtains 

5(53 1105^6 

^^ = ^°+ 8^-12^^ + ^^^^' (^) 

where 5^ = [m^'Y/A. It is known that this constitutes a semiconvergent asymptotic series 

solution to (^. 

2 Semiclassical propagator 

In order to study the effect of a space-time varying mass m^ on the propagator, we use the 
Schwinger-Dyson equations in the Keldysh closed time contour (CTC) formalism [0, H. This 
implies the following formally exact propagator equation in the Wigner representation 

cosO{fc2-m2±icjr}{G"''^} = l. (5) 



where for simplicity we neglected the self-energy contribution. The O-operator denotes a 
generalized Poisson bracket 0{f}{g} = {d^fdkg — dkfdxg)/2 and fc^ = {u, k) is the canonical 
momentum and x^ = {t, x) is from now on the average coordinate. The propagator equa- 
tion in the Wigner representation (|^) is useful when the background is slowly varying, i.e. 
formally when d^OkG is much smaller than G. In our notation C"'" indicates the retarded 
and advanced propagators, whose damping terms are izdjF, respectively. To the 0th order 
in gradients the solution to Eq. (^ reads 

GT = 17, ^^-p ■ (6) 

From Eq. (^ we then infer that the gradient correction occurs first at the second order P]. 
The propagator equation @ defines the spectrum of excitations in the system. This 
information is encoded in the spectral function A = ^ (G"" — G") which, as a consequence of 
the general equal time commutation relations, satisfies the well known sum rule: 

1 r°° 

- / dLULuA = 1 . (7) 

n j-oo 

Using the complex properties of the lowest order propagator (|^), one can rewrite the sum 
rule as a contour integral along the path shown in Fig. 1 of the complex valued propagator 
Go = (^^ — ^o)^^' "where u^ = k"^ + m? and we set F — > 0. It is then easy to see that 
(i/vr) f(.^ rfcucjGo = X]w=±wo -^^^P^^o] = 1- It is convenient to introduce a new variable 
z = k"^ — m^ . In terms of z the sum rule for Go = I/2; becomes simply 

-^ / dzGo = Res[Go].=o = 1 , (8) 

zvr JCz 

because one can deform the contour C^ to C^ around the origin, as shown in Fig. 2. Note that 
the mapping oj ^> z = u"^ — u"^ induces a double covering of the complex 2;-plane; one for 
particles (Re [oj] > 0), and one for antiparticles (Re [uj] < 0); Re [u] = is mapped onto the 
negative real 2;-axis. Strictly speaking the two covers are identical only when the propagator 
conserves CP symmetry. 

The sum rule can also be computed by integrating along a circle C/j at R^ 00, shown 
also in figure 2. We now assume that this is true in general for the full physical propagator 
G. This means that G must have a residue equal to unity at |2;| =00, which is equivalent to 
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Figure 1: The integration contour for spectral integral in complex w-plane. 
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Figure 2: The integration contours for spectral integral in the complex z-plane. 



the following simple boundary condition 

z 



when \z\ -^ oo. 



(9) 



More generally we will assume that any spectral integral of the form / dioLoAT ^ where T is 
some dynamical quantity, is to be computed by a contour integral of GT over a contour C 
which contains contributions from all poles and possible discontinuities of the propagator. 

As was already pointed out, the first nontrivial correction to the propagator equation (^) 
occurs at the second order in gradients. To this order one finds the following equation for 
the complex valued G: 



{e-rr?)G^\[dl^'^rr?"dl 



G 



(10) 



where for simplicity we considered a stationary case vr? = 'm^(x), again with a;^ = (0, x, 0, 0). 
Assuming that the derivative corrections are small, it is appropriate to use an iterative 
procedure, which to the lowest nontrivial order gives 



G, 



6, 



5i 



;ii) 



where 82 = vn}" jl^ and 6\ = {m?'Y /2 + kim?" . A similar iterative procedure was used 
in Eq. (^), when constructing the WKB approximation to the wave function. In the rest 
of this letter we explain in what sense this iterative procedure can be interpreted as the 
semiclassical approximation to the propagator. 

Let us mention that the multiple pole form of the propagator ([TID can, in a loose sense, 
be reconciled with the picture of quasiparticle excitations. To see this, note that to the 
second order in gradients G2 can equivalently be written as 

G'2= ," ^, ^E^^. (12) 

where a„ denote the residues at the simple poles Zn- These residues and poles are in general 
complex however, yielding unphysical complex dispersion relations z = Zn and weights a„. 
Indeed, the true physical meaning of the semiclassical propagator emerges only when G is 
understood in an operational sense inside the spectral integral over the complex path in the 
sense conjectured above. 

3 Recovering the WKB result 

To study the physical consequences of Eq. (j^) in more detail, without overly complicating 
the calculations, we now make the same simplifying assumption that was used to obtain the 
equation (Q): d'^m? = for n >2. Without further approximations Eq. (j^) then becomes 

6'^ + zG = l, (13) 

where Si/2 ^ 6^ = (m'^'Y/A and z = k"^ — -m^. The standard solution of (|13D is known. 
The associated Airy function Hi(z) ||^ is unphysical however, because it does not satisfy the 
boundary condition (^ and it has no isolated poles, making it impossible to satisfy the sum 
rule. However, we can obtain a formal series solution with the desired properties through an 
iterative procedure: 

N^ooz^^^ ' 3"n! \zj ^ ' 

This can further be formally resummed to yield 

G = lim - / rfMe-"-^^"'^/")' . (15) 

A^oo z JO 



In the limit A ^ oo, G is in fact equivalent to the associated Airy function for positive 
real z, but it diverges for Re [z^] < 0. This divergent behaviour is to be expected due 
to the unbounded nature of Eq. (0), and it is also present in the series solution, which 
is formally divergent for any finite \z\. However, we show below that taking a finite A in 
(p!5|), or truncating the asymptotic series in Eq. (|l^ to some finite iV, leads to a meaningful 
regularization, and in fact to the best estimator for the propagator. 
To this end consider spectral integrals of the form 

2 /""^ ., ,_ if dz 

M 



G[T] ^ - / dKAT ^ — i -^^GT, (16) 



^0 



where z = k^ — k^ with k^ = uo — k^f — m , and the contour C^ is the one shown in figure 
2. The test function T representing a dynamical quantity can be any meromorphic function 
with no poles on the real tu-axis. Note that, as a consequence of separability of spectral and 
dynamical properties, this is true for the important example of the generalized distribution 
function in the dynamical Schwinger-Dyson equation. Furthermore, T must admit a Taylor 
expansion around the propagator poles, i.e. it admits gradient expansion. 

We chose integration over momenta rather than frequencies in Eq. (|T6|) to facilitate direct 
comparison with the standard form of the semiclassical result in Eq. (^. Expanding T into 
the Taylor series in k around k = ko, and computing the residue one obtains 



G[T] = Res 



GT/Jk, 



= T.tnrt^ (17) 

2=0 n = 



where Tq = T{ko), Tq = {dk^T)k^=koi etc. It suffices to consider the first two terms in (|T 
to obtain the semiclassical dispersion relation. Indeed, 

G[r] = ^ + (i-f]n + o{r") 



k 
where 



^(^sc) ^ ^^^„^ 



.._ig^_^,„rp^^/.y; ^^^, 



It can now be easily verified by inspection, that fcgc is equivalent to the WKB momentum 
kx in Eq. (1). Moreover, one finds that under the assumption d'^m? = for ri > 2, Eq. (0) 

6 



reduces to the following differential equation: y'" + C,y' + y/2 = for y = S^^'^/kgc, where 
( = kl/6 and y' = dy/d(, also satisfied by the solution ([T9|) . 

In figure 3 we plot the relative deviation k, = ksc/ko — 1 from the unperturbed momentum 
kf) vs. S/k^ for the three lowest order solutions (2nd, 4th and 6th order in gradients) as well 
as the formally resummed integral representation 

hl= r ^ e-^-\iu8/klf _ (20) 

k^r Jo \/t^ 



For small S/k^ the agreement is extremely good and it becomes progressively worse as 
S/kQ approaches 0.3. It can be shown that the highest possible accuracy of the order 
exp — (2/3)(A;q/5)^'^ is reached by the propagator G2N truncated at the 2A^-th order in deriva- 
tives, where A^ is defined by 5 ~ A2Ar = kQ{l/3Ny^^. In fact the integral expression is just 
the envelope of G'2Af's at 5 ~ A2Ar, and hence represents the best possible estimate of fcsc 
for any value 5 < A2. However, even though the integral gives a reasonable result also for 
large values (5 > A2, one should be careful in attributing a physical meaning to it, because 
analytic continuation of an asymptotic series is not unique. 

It is interesting to note that in the case of the electroweak phase transition the expanding 
phase transition front typically has a width ~ 10 — 2O/T0, so that typically S/k^ <, 0.1. In 
this region the asymptotic series is accurate to 10 decimal points. 

We have seen that the semiclassical spectral function projects a dynamical quantity T 
"on-shell" by shifting the energy hypersurface from k^ to ksc = ksc{uj,kf,x). This projection 
however induces additional structure in Eq. (|T^ contained in the higher order derivative 
terms 7^", etc., and hence the semiclassical limit in field theory cannot be described by a 
single on-shell function and the semiclassical dispersion relation. This is a direct consequence 
of the higher order poles in the semiclassical propagator, which occur generically in the 
presence of a space-time varying background. Note that only a finite number of terms in the 
expansion in 7^ in Eq. ([T6| ) contribute. Indeed, the coefficient of the n-th order derivative 



Tq /n\ has the form {JkQ — z — kQ)'^ = {—l)^{z/2kQ)^ + 0{z'^^'^), and hence the asymptotic 
series representation (jM]) of G implies that, at the order 2N in gradients, only the elements 
Tq with n < 3N can have nonvanishing coefficients. For example, using the second order 
approximation to the propagator in Eqs. (|lT]-|l^), the series in Eq. (^) has only 4 terms. 




Figure 3: Plotted is the quantity n = ksc/ko — 1 for three lowest order gradient approxima- 
tions for G along with the formal integral solution (solid line). Inset: deviation of /€„ from 
the integral result for small S/k^. 

terminating with ((5'^/24A;q)7^"'. The meaning of these corrections can be fully appreciated 
only when studying the semiclassical approximation of concrete dynamical quantities. Of 
particular importance is the quantum Boltzmann equation, which can be obtained from the 
dynamical Schwinger-Dyson equation by the method of spectral projection described in this 
letter. In this case the test function is related to the density matrix of the system and so 
the higher order derivative terms represent nonvanishing off-diagonal correlations between 
different momentum states. This is just a reflection of the uncertainty principle, according 
to which localisation in the position necessarily implies delocalization in the momentum. 

4 Conclusions 



In this letter we studied the semiclassical limit in field theory. We found that in the pres- 
ence of slowly varying backgrounds, there is a simple iterative procedure by which one can 
construct the semiclassical propagator as a semi-convergent asymptotic series controlled by 
powers of gradients. While we discussed at length the simple case dl^m'^ = for n > 2, 
the method can straightforwardly be extended to the general case with higher derivatives 



and nonvanishing self-energy. The semiclassical dispersion relation appears in a nontrivial 
way through study of spectral integrals and cannot be associated to a simple shift of the 
pole of the propagator. We further found that the spectral projection of a dynamical quan- 
tity necessarily contains higher order momentum derivatives as a reflection of localization 
in coordinate space. As a consequence the plasma dynamics is not described by a single 
classical Boltzmann equation, but instead by a set of independent equations describing not 
only evolution of the distribution function, but also of a limited number of its derivatives. 
We will pursue this issue in a forthcoming publication. 
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